674 



Ind. Eng. Chem. Res. 1993, 32, 674-684 



Gas Absorption Studies in Microporous Hollow Fiber Membrane 
Modules 
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Department of Chemistry and Chemical Engineering, Stevens Institute of Technology, Castle Point, 
Hoboken, New Jersey 07030 

A comprehensive experimental investigation of gas-liquid absorption in a ahell-and-tube type 
microporous hydrophobic hollow fiber device in a parallel flow configuration was carried out. Two 
modes of countercurrent gas-liquid contacting were studied, the wetted mode (absorbent liquid 
filled pores) and the nonwetted mode (gas-filled pores). The absorbent flowed through the fiber 
bore in most of the experiments. The systems studied include pure C0 2> pure SO2, COr~Nz mixtures 
and SOa-air mixtures. The absorbent was pure water. The absorption process was simulated for 
each case with a numerical model for species transport with and without chemical reaction. Laminar 
parabolic velocity profile was used for the tube-side flow, and Happel's free surface model was used 
to characterize the shell-side flow. The model simulations agreed well with the experimental 
observations in most cases. SO2 removals as high as 99 % were obtained in small compact contactors. 
High Kua and low height of transfer unit (HTU) values were obtained with hollow fiber contactors 
when compared to those of conventional contactors. 



Introduction 

In conventional large-scale gas absorption processes, the 
gas mixture is contacted with a scrubbing liquor in devices 
like packed towers, spray towers, venturi scrubbers, and 
bubble columns. These dispersion-based contacting meth- 
ods have many drawbacks. An approach having much 
promise and without these shortcomings is membrane- 
based contacting using microporous hydrophobic hollow 
fibers. The microporous membrane-based device acts as 
a gas absorber, with the gas flowing on one side and the 
absorbent liquid flowing on the other side of the membrane 
without either phase being dispersed in the other. The 
gas-liquid interface is at the pore mouth. Unlike that in 
a true liquid membrane process, the absorbent liquid has 
to be regenerated in a separate stripper which can be also 
a membrane-based contactor. 

A membrane-based contacting device can have a number 
of configurations. The focus here will be primarily on 
parallel flow modules of microporous hollow fibers. The 
membranes can be hydrophobic or hydrophilic, and the 
pores can be either gas filled or absorbent filled. The 
absorbent may flow on the tube side or shell side. These 
membrane-based gas separation devices offer several 
advantages over conventional contacting devices: much 
higher surface area per unit contactor volume; independent 
control of gas and liquid flow rates without any flooding, 
loading, weeping, or foaming; known gas-liquid inter facial 
area; modularity, easy to scale up or down; availability of 
all of the membrane surface area for contacting regardless 
of how low the individual phase flow rates are. 

These advantages have led a number of investigators to 
rplore different applications. The first known application 
of a microporous membrane as a gas-liquid contacting 
device using hydrophobic flat Gore-Tex membranes of 
polytetrafluoroethylene was for oxygenation of blood 
(Esato and Eiseman, 1975). Tsuji et aL (1981) used 
hydrophobic microporous hollow fibers of polypropylene 
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Figure 1. Nowetted mode of hydrophobic membrane-based gat- 
liquid contacting. 

for blood oxygenation. The applications of direct interest 
here are those for acid gas cleanup. Qi and Cussler 
(1986a,b) studied the absorption of a number of gases using 
various solvents and aqueous solutions in a hydrophobic 
microporous hollow fiber device where the membrane pores 
were gas filled. Yang and Cussler (1986) have noted a 
10-fold increase in the gas separation efficiency when using 
a hollow fiber device over a conventional packed column. 
Ogundiran et a! (1988, 1989) in their experiments with 
absorption of SO2 using microporous hydrophobic hollow 
fibers found this technique to be an extremely promising 
alternative to conventional scrubbers used for flue gas 
desulfurization. Sirkar (1992) has recently reviewed 
microporous membrane-based gas absorption processes. 

The above-mentioned contacting studies were carried 
out using a microporous hydrophobic membrane. Figure 
1 shows an absorbing liquid (that does not wet the fibers) 
flowing on one side and gas flowing on the other side. The 
gas pressure (P go ) has to be lower than that of the liquid 
(Paq) to prevent dispersion of gas as bubbles in the liquid. 
Liquid does not enter the pores unless a certain critical 
pressur is exceeded (Kim and Harriott, 1987). This mod 
of operation is identified as the nonwetted mode since the 
liquid does not wet the membrane and the membrane pores 
are gas filled. The gas/liquid interface shown in Figure 
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Figure 2. Wetted mode of membrane-based gas-liquid contacting. 

1 lies at the pore mouth of the hydrophobic fiber on the 
absorbing liquid side. In this mode the gas species being 
separated diffuses first through the gas in the membrane 
pores before dissolving in the absorbent. Because the pores 
in the fibers are gas filled, the gas-phase mass-transfer 
coefficients are much higher than the liquid-side mass- 
transfer coefficients if no fast reaction is taking place or 
when the gas solubility is very low. Thus gas absorption 
is likely to occur without practically any increase in mass- 
transfer resistance due to the presence of the membrane. 

There is yet another mode of operation where the 
membrane pores are absorbent liquid filled; this is called 
the wetted mode and is shown in Figure 2. Here a 
hydrophilic membrane is spontaneously wetted by the 
absorbing liquid which flows on one side, while the gas 
flows on the other side. In the case shown in Figure 2, the 
gas pressure has to be higher than the liquid pressure to 
prevent the liquid from dispersing as drops in the gas. A 
hydrophobic membrane can also be used in this mode using 
a nonwetting liquid by employing an exchange process 
(Bhave and Sirkar, 1986, 1987) which allows one to 
incorporate a nonwetting absorbing solution in the hy- 
drophobic membrane pores. This mode may be useful 
when there is a very fast or instantaneous liquid-phase 
reaction and the gas resistance would be controlling the 
rate of transfer. The membrane-phase resistance may be 
eliminated in such a mode of operation. There is another 
reason for investigating such a mode. If long-term 
operation of a hydrophobic membrane leads to wetting of 
the pores, such a mode of operation may be adopted. 

In this paper, the absorption of acid gases like S0 2 and 
CO2 have been investigated using modules of microporous 
hydrophobic hollow fine fibers. Experiments have been 
carried out with both gas-filled membrane pores (non- 
wetted mode) and liquid-filled membrane pores (wetted 
mode) using pure gases and their mixtures, with different 
absorbing liquids flowing either on the tube side or on the 
shell side. This paper focuses on pure water as absorbent 
A second paper will cover diethanolamine in water. 
Although the absorption of C0 2 in water is essentially 
physical, absorption of SO2 in water is very strongly 
influenced by ionization so that there is reactive absorption. 
Gas absorption models have been developed on the basis 
f first principles to explain the observed separation and 
evaluate some commonly used correlations. The present 
analysis attempts to pr vide a fundamental basis for 
understanding aspects of this promising contacting tech- 
nique. The experimental values of K\ja and HTU btained 
here have been compared with those available in the 
literature. Gaseous systems investigated here include pure 
CO2, pure SO2, CO2-N2 mixtures, and SOr-air mixtures. 



Experimental Section 

Micropor us hydroph bic Celgard X-20 and X-10 
polypropylene h Uow fibers made by Hoechst Celanese 
Corporation (SPD, Charlotte, NC) were used. Able to 
withstand air temperatures up to 75-80 °C, these fibers 
have an internal burst strength greater than 220 psi. Both 
fibers have a bubble p int pressure of 150 psi and a mean 
pore diameter f 0.03 urn- The characteristics of hollow 
fiber modules 1-4 fabricated in the laboratory and used 
here are provided in Table I. Module 2 has the highest 
area/volume. 

The gas composition was measured by a Hewlett* 
Packard 5890A gas chromatograph using a thermal con- 
ductivity detector (TCD). The detector response was 
recorded by a digital reporting integrator 3393A (H J 5 ., 
Paramus, NJ). A 10-port gas sampling valve was utilized 
(Sengupta et aL, 1990). A two-column system based on 
a S-ft X Vr in. Chromosorb 108 column (Chrompack Inc., 
Bridgewater, NJ) and a 10-ft X 1 /$-in. molecular siev 
(All tech Associates, Deerfield, IL) column was used. The 
column sequence was reversed after 1.9 min to separate 
O2 and N2 by the molecular sieve column. By this reversal 
CO2, SO2, and moisture did not enter the molecular sieve 
column thereby ensuring it a longer lifetime. Sampling 
was done continuously. 

The liquid-phase concentrations were analyzed and used 
to compare the gas analysis via mass balances. Inthecas 
of pure CO2 if the change in gas flow rate was limited, 
concentrations based on the liquid phase were considered 
more reliable. For CO2, liquid-phase-concentration mea- 
surements were made using a CO2 electrode (Orion Model 
94-02, Orion Research, Cambridge, MA). The measure- 
ment range of this electrode was 4.4-440 ppm CO2 at a pH 
range of 4.8-5.2. Reproducibility was ±2%. Titration 
procedures were also used to cross-check the data (Betz, 
1962). The concentrations measured by the electrode and 
the titration method were in agreement in the low rang 
0-100 ppm. In higher concentration ranges, a different 
titration method was used (Snell and Hitton, 1966). In 
the case of SO2, the method of analysis was based on the 
procedure described in Standard Methods for the Ex- 
amination of Water and Waste Water (AWWA, 1965). 

Figure 3 shows the experimental setup used for gas 
absorption studies. The gas mixture to be separated or 
the pure gas to be absorbed was passed through the modul 
either on the shell side or on the tube side, and th 
absorbing liquid flowed countercurrently through the 
absorber on the other side, separated from the gas by the 
membrane walls. The gas flow rate was controlled by an 
electronic mass flow transducer and a flow controller 
(Models 8100 and 8200; Matheson, East Rutherford, NJ). 
Pressure gauges at the module inlet and outlet indicated 
the gas pressures which varied between 0 and 14 psig 
depending on the mode of operation. The gas stream outlet 
concentration was measured by a gas chromatograph (GC), 
after which the gas was vented into the laboratory exhaust 
hood. A needle valve was used to control the liquid flow 
rate. The liquid flow rate was measured at the outlet of 
the membrane contactor at regular intervals. The back- 
pressure regulator (BPR) was not needed when the 
membrane pores were gas filled; i.e., the system was 
operated in the n nwetted mod . Experimentally ob- 
served concentration change from the feed to the outlet 

f the module in the gas phase was used to determine the 

verall logarithmic mean based verall mass-transfer 
coefficient Koim, using the f Uowing equation: 
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Table I. Ch*racterif tice of Hollow Fiber Module* Ueed In This Study 



module 
no. 


fiber o.d. t 
Mtn 


fiber i.d., 
jim 


shell diam, 
cm 


length of 
module, cm 


no. of 
fiben 


porosity 


packing 
fraction 


area of 
contact i.d., cm 2 


o, 

cm 2 / cm 3 


1 


150 


100 


0.4672 


30.0 


102 


0.2 


0.1098 


96.13 


19.52 


2 


300 


240 


0.4927 


30.0 


132 


0.4 


0.4887 


298.57 


52.2 


3 


300 


240 


0.4672 


18.0 


75 


0.4 


0.3229 


101.78 


34.44 


4 


300 


240 


0.4672 


17.6 


76 


0.4 


0.1098 


98.96 


34.44 




T GAS OUTLCT 



Figure3. Gas absorption using hollow fibers: experimental setup. 



01m SQA^Cte 



(1) 



Kohn is also denoted as Kt in this paper. Ay is based on 
the surface area of gas-liquid contact The quantity AG m , 
the logarithmic mean concentration difference of gas 
species, is given by 



AC, 



lm 



(HC g in - C\ out ) - (HC g out ) 

H Hc gout ; 



(2) 



Here H is Henry's law constant which is the ratio of the 
interfacial concentration of absorbed species i in liquid 
divided by that in the gas: 



*HCu 



(3) 



if species i is represented by A and the liquid phase by w 
for water. In the case of CO2-H2O system, the value of 
Henry's law constant is 0.8314 at 25 °C (Geankoplis, 1972) 
and that for S0 2 -H 2 0 is 26.86 at 25 °C (Hikita et el, 
1978). 

Model Development 

A theoretical model is developed here to describe the 
absorption of pure CO2, pure SO2, and their mixtures in 
nitrogen or air by water. Both wetted and nonwetted 
modes of operation are considered here, and in both modes, 
the liquid flow was always countercurrent to the gas flow 
and on the tube side. The following assumptions were 
utilized: (1) steady state and isothermal condition; (2) no 
axial diffusion; (3) fully developed parabolic liquid velocity 
profile in the tube side; (4) use of Happel's free surface 
model (1959) to characterize the shell-size velocity profile; 
(5) ideal gas behavior; (6) applicability of Henry's law. 

Using these assumptions, the general f rm of the 
conservation equation for the absorption processes is given 



nuc SURFACE 




FIBER WALL / 



2 TUBE-SIDE " ) 
LIQUID FLOW 1 




JHELLjSIDE CAS _ 






2=0 2=L 

Figure 4. Free surface model (adapted from Gill and Banna! ( 1973) ). 



(4) 



In the above equation, i denotes the components involved 
(i « A, B, .„), ; denotes the phase (gas (g), liquid (1)), and 
k denotes the location (shell (s), tube (t), or membrane 
(m)). The quantity v t k is the prescribed axial velocity 
profile for the tube or the shell side. For the tube side a 
laminar parabolic velocity profile is assumed: 

u lt = 2< y „>[l-(r/r t ) 2 ] (5) 

The flow in the shell side is usually very complex. Happel's 
(1959) free surface model as represented in Figure 4 is 
used to characterize the shell side velocity profile; it is 
given by the following equation: 



u„(r) » 2< 



r (r/rf-(r 0 /r/ + 2]n(r 0 /r) 1 
L 3 + (r^rf - 4(Vr e ) 2 + 4 lnfr^r,) J * 



where r # is the radius of the free surface and is given as 

[1 "h/2 

Pure Gas Absorption: Nonwetted mode. In this 
mode of absorption of pure C0 2f the pores of the membrane 
are tilled with pure gas flowing in the shell side. There 
is no gas-phase-transport resistance. Due to the low value 
of the equilibrium constant (Ke ■ 4.4 X 10~ 7 mol/L), 
bicarbonate formation is neglected. Equations 4 and 5 
are solved for the liquid phase for i ■ A ■ C02, j m l t k 
« t, and r A = 0, along with the following boundary 
conditions. The concentration profile is symmetrical about 
the center of the tube and therefore satisfies 

0 (7) 

Further, since pure CO2 gas of known pressure is present 
in the shell side, the species concentration at the inside 
wall (r * r t ), is known. Hence the boundary condition at 
the wail is 

C ih (r t S) - C iw (8) 

Here C IW is the product f Henry's law constant times 
gas-phase concentration and is given by eq 3. In most of 
the experiments conducted, the pure gas pressure remained 



essentially constant along the module length. Therefore 
C,w remained constant with z (which is essentially the 
Graetz solution (Skelland, 1974)). There were a few 
experimental situations where the pure gas pressure 
changed somewhat along the module. In such a case i C; w 
varied with axial location z (linear variation was assumed 
for numerical calculations). The inlet condition is also 
known 

C m (r,0)«C inlit «0 (9) 

So far in the species conservation equation, there was 
no reaction. The situation changes for absorption of pure 
SO2 in water. The following hydrolysis reactions take place 
when SO2 is absorbed in pure water: 

S0 2 + 2H 2 O^H 3 0* + HS0 3 -, K E = 0.014 M (10) 
HS(V + H 2 0 p* H 3 0 + + S0 3 2 -, JC E - 6.24 X 10"* M 

If the pH is less than 4-5, formation of S0 3 2_ is negligible. 
By choosing operation conditions that exclude the pH 
range mentioned above, S0 3 2 " formation is avoided 
(Seinfeld, 1980; Roberts, 1979). In the present study the 
operational conditions were such as to allow the assumption 
that only HSO3* formation is effective. 

The diffusion coefficients of HS0 3 - (species B) and H 3 0 + 
at 25 °C are 1.47 x 1(H> and 9.3 X 10- 5 cm 2 /s, respectively. 
An electrical potential develops which retards the motion 
of HaO+ and speeds up HSO3-, so as to maintain electro- 
neutrality. Thus an effective diffusion coefficient in the 
liquid is defined, which in the binary case is given as 
(Vinogard and McBain, 1941) 



'Bl ' 



^H ) O^HSO a - 
^H,0* + ^HSO,- 



(12) 



The rate expressions for the reaction (eq 10), after using 
the electroneutrality condition, Chso," ■ Ch,o*i are ex- 
pressed by the following equations: 



■ r so, : 



so 8 



+ ^(CHSO r ) 2 



(13) 



(14) 



where the forward rate constant k x is given as 3.17 X IQr 2 
s' 1 (Wang and Himmelblau, 1964). Since eq 4 is valid for 
both components SO2 (A), and HSO3" (B), there are now 
two coupled partial differential equations to be solved 
along with the boundary conditions. With eqs 13 and 14 
representing the reaction terms -fx and -rs, respectively, 
eq 4 for SO2 and eq 4 for HSO3- are solved simultaneously 
along with the above boundary conditions (eqs 7 and 8). 
In the case of HSO3- at the wall, equilibrium is assumed 
and the following equation is used to represent the wall 
condition: 



1/2 



(15) 



For details, see the thesis by Karoor (1992). 

Wetted Mode. In this mode of operation the liquid is 
flowing in the tube and is also incorporated in the 
membrane pores. An additional equation is required to 
describe the diffusion through the membrane pores as the 
liquid-filled pores offer additional resistance. At high 
porosity (greater than 10% ), Keller and Stein (1967) have 
shown that the diffusional process is essentially one 
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dimensional. The following assumptions were made in 
deriving the membrane-phase equation: 

1. The effective membrane-phase diffusion coefficient 
Am i8 defined by Dat/r. 

2. The diameter of the pore is uniform. 

3. There are no gradients in the direction normal to the 
pore axis. 

4. There is no convective transport. 

The species balance eq 4 is still valid for the liquid in 
the membrane phase (j : « 1, k - m) for each species without 
the convection term. The following boundary conditions 
are used at the membrane-tube interface: 



(16) 



A flux continuity at the same boundary is assumed, giving 
the equation 

The equation for the membrane phase is solved using the 
above boundary conditions, namely eqs 16 and 17, to yield 
the following equation at the inner wall of the membrane 
r « r t (Karoor, 1992). 

W^-C^f Qnr^tor^ (18) 

Now eq 4 for i = A, / = 1, and k ■ t along with eq 18, the 
boundary condition eq 7, and the wall condition at r « r 0 , 
Cmn(ro^) = C,w are solved. 

The above system of equations is for CO2, since there 
is no reaction in the liquid phase. For S0 2 , two species 
SO2 (A) and HSO3" (B) are present and eq 4 is valid for 
each species for the tube side. The equations for the 
membrane phase are given as 



Ala j£ U Blm 



A flux continuity equation at the boundary is assumed, 
giving the equations 



r»r„ 



r-r u 



Blm, 



831 dr 



(22) 



Equations 19 and 20 are solved along with the boundary 
conditions eqs 21 and 22 and then combined with the 
equations for the interior and the boundary conditions. 

Absorption from a Gas Mixture. In gas absorption 
from a mixture, the gas-phase resistance has to be 
considered; an additional equation for the gas phase has 
to be solved. For CO2 absorption in water, eq 4 is also 
valid for the gas flowing in the shell side. The shell-side 
gas velocity profile is given by eq 6. At the free surface 
r -> r«, the following boundary condition is valid: 



(23) 



At two membrane interfaces the following conditions are 
valid: 



Jr." H 



(24) 
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Table II. Parameter* Used In the Absorption of SO* and CO* in Water 



temp 
•C 



solubility 0 of 
mol/liter 



diffusivity 0 
of S0 2 d&D), 
cmVs 



solubility* 
ofCOa(tfm). 
mol/L 



diffusivity* 
of CO* (l&D), 
cmVi 



equilib* constant 
ofSOj (K), 
mol/L 



16 

25 



1.5420 
1.0559 



1.34 
1.76 



0.045446 
0.3398 



1.40 
1.92 



2.19 
1.70 



° Hikita et al (1978). * Geankoplii (1972). « Majumdar et aL (1988). 

' = C Ue -C U t {25) 

Flux continuity equations at the same boundaries are also 
valid, giving the equations 



r = r tl 



r = r 0l 



(26) 



(27) 



Then, eq 4 is solved for the membrane phase with the 
above boundary conditions to give the gradients at r ■ r t 
and r = ro, as 



*S _ \ . C MDk» 1 



[In r 0 - In r t ]r 0 



(28) 



(29) 



The equations for the shell and tube are combined with 
the above equations for the membrane inner and outer 
boundaries and solved simultaneously with the boundary 
conditions eqs 7 and 23 and the inlet conditions for CO2 
absorption. In the case of SO2 absorption all the above 
equations and boundary conditions are valid except that 
the reaction term has to be taken into account; in addition 
similar relations are written for HSO3-. At the gas-liquid 
interface the formation of HS0 3 ~ is assumed to be zero 
(for details, see Karoor (1992)). 

The above equations are solved using a finite difference 
scheme (Davis, 1984). The diffusion terms are discretized 
using the centered in space scheme, whereas the convection 
term is discretized using the upwind difference scheme. 
A downward difference scheme is used when the velocity 
in the shell side is in the opposite direction to that of the 
tube side. Using the method of lines, the above partial 
differential equations were reduced to a set of algebraic 
equations. These equations are tridiagonal in form as a 
result of the discretization method adopted. This set can 
now be easily solved using the Thomas algorithm (Karoor, 
1992). In the case of SO2 absorption the equation is solved 
for Cxit (based on Ceit values determined in the earlier 
iteration). The equations for Con are solved using a similar 
approach; the nonlinear term in eqs 19 and 20 are 
pseudolinearized by representing Ceim 2 and CbizdCbud*! 
where Caim* is the value of Cbi» from the previous iteration. 
The solution procedure is continued till convergence is 
attained. 

Results and Discussion 

The results of countercurrent absorption of pure CO2 
and pure S0 2 in water flowing in the tube side of the hollow 
fiber modules are presented and discussed first. Perfor- 
mance of gas mixtures is provided later. Th values of the 
parameters used are given in Table II. The experimental 
results have been compared with the results btained from 
the numerical solution of the governing partial differential 
equations. Limited p rformance data have been provided 
also for shell-side water flow. Both experimental and 
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Flfure 5. Comparison between wetted and nonwetted modes of 
operation for absorption of pure CO] in water. 

model results have also been recast in the fashion of a plot 
of Sherwood number vsGraetz number. Besides the pure 
gas absorption experiments, the experimental results and 
the model simulations are discussed for cases where the 
feed gas is either a mixture of C0 2 (40%, 10%) in N 2 or 
a mixture of SO* (1%, 0.95%, 0.4%) in air. One or two 
or all three of the hollow fiber modules 2, 3, and 4 were 
used in these studies. 

Figure 5 illustrates the variation of the liquid outlet 
concentration of CO2 as a function of water flow rate in 
module 2. The gas flow rate was held constant at 76 cm 3 / 
min for all the runs shown in the figure as the liquid flow 
rate was varied. Data and numerical simulation results 
are presented for both nonwetted mode and wetted mode. 
It is clear that the liquid outlet CO2 concentration in the 
wetted mode is lower than that in the nonwetted mode. 
This effect is more apparent as the liquid velocity is 
increased. Since the liquid boundary layer resistance 
decreases with increasing liquid velocity, the membrane 
resistance starts to control. By filling the pores with the 
absorbent liquid, the membrane-phase resistance is in- 
creased as the liquid offers additional resistance to the 
transport; this behavior highlights the fact that this is a 
liquid-phase-controlled process. Hence it would not be 
advantageous to operate under wetted-mode conditions 
for this particular system. Exactly the same trends were 
observed in the case of module 3. The theoretical curves 
shown in Figure 5 were obtained by solving eq 4 along 
with the boundary conditions. The theoretical predictions 
for these situations have a maximum deviation of 7 % from 
the experimental results. 

The effect of using a module with a higher contact area 
can be observed by using a greater number of fibers or 
increasing the length of the module. Figure 6 shows a plot 
of the concentration f CO2 in the liquid outlet f r three 
different modules of varying surface area as a function of 
the liquid flow rate in the nonwetted m de. From Table 
I it is seen that module 2 has the maximum c ntact area 
followed by module 3 and then module 4. In the case of 
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Figure 6. Comparison of concentration of C0 2 in the liquid outlet 
for modules of different surface area. 
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Figure 7. Sherwood number versus Graetz number variations for 
pure COs absorption in water with different modules. 

module 2, a higher amount of CO2 in the liquid outlet is 
observed for the same liquid flow rate when compared 
with modules 3 and 4. Obviously using a longer module 
would require a higher pressure drop to achieve a given 
flow velocity. A maximum CO2 removal of 76% was 
achieved in module 2 with a liquid flow rate of 244 cm 3 / 
min. The numerical model simulations describe these 
behavior quite well 

Figure 7 illustrates the same experimental results for 
modules 2, 3, and 4, but now as a plot of Sherwood number 
versus Graetz number, where Sherwood number is defined 
as Koin>(2rt)/D A t f° r ^ e nonwetted mode (for the wetted 
mode r t is replaced by r 0 ) and Graetz number for the tube 
side is defined as Nr»Nsc rrJ2L. Ko\m is the logarithmic 
mean overall mass-transfer coefficient calculated by eq 1. 
A logarithmic mean driving force is used in the calculations 
as given by eq 2. An increase in Graetz number due to an 
increase in the the liquid flow rate increases the mass- 
transfer coefficient as expected. Another feature that can 
be observed is the effect of operating the module in the 
wetted mode. In the wetted mode the pores of the fiber 
are tilled with th abs rbent liquid. This increases the 
resistance to mass transfer across the membrane p res, as 
a result of decreased gas diffusi n coefficient in the liquid 
phase. This effect is clearly sh wn in Figure 7 via a 
decrease in the Sherwood number by almost an rder f 
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Figure 8. Effect of tortuosity on the CO2 outlet concentration for 
module 2. 

magnitude. As is evident from Figure 7, the effect of the 
wetted mode in module 3 was also observed to be similar. 
The extent of decrease in mass transfer in the wetted mode 
in module 3 is also comparable to that of module 2. 

Figure 7 also illustrates the results obtained from the 
theoretical model for the nonwetted mode. The model 
results and the experimental observations show good 
agreement with a maximum deviation of about 11.5%. 
The Leveque solution (Leveque, 1928), also plotted in this 
figure, shows a reasonable agreement with the data. The 
Leveque solution results are somewhat higher at higher 
Graetz numbers. 

In the wetted mode, an effective membrane diffusion 
coefficient was defined as 2) Am = £W/t, where c/r was 
taken to be 0.4/2.5, from the experimentally measured 
values by Prasad and Sirkar (1988) for X-20 fibers with 
240-Mm inner diameter and 300-/*m outer diameter. In 
the case of X-10 hollow fibers with 100-Mm inside diameter 
an 150- /an outside diameter, the quantity t/r was taken 
to be 0.2/3.0. 

Using this definition of effective diffusion coefficient 
and the values of e/r, the model discussed before for the 
wetted mode was solved numerically to obtain the results 
shown in Figure 7. The agreement between the model 
and experimental results is good with a maximum deviation 
of 7.5 % for module 2. For module 3, Figure 7 shows that 
the agreement between the model results and experiments 
is again good with a maximum deviation of 12%. As 
discussed before, the length of module 2 is higher than 
that of module 3; also the number of fibers is greater in 
module 2 when compared to module 3. Thus, under similar 
operating flow rates, the Graetz number is lower for module 
2 when compared to module 3, and it results in lower 
Sherwood numbers for module 2. 

In Figure 8 the predicted CO2 outlet concentration for 
module 2 in the wetted mode is plotted by allowing a 
maximum variation of ±20% in the tortuosity r for X-20 
fibers (2.5) used in the present study. As the tortuosity 
increases, the transfer of CO2 into the liquid decreases 
since CO2 has to diffuse through longer lengths of liquid. 
If the measured tortu sity is ±20% different from the 
value 2.5 used, then the maximum difference in the 
concentration is around 18 % for the liquid flow rates and 
the system under study. 

Figure 9 illustrates the results for module 3 when pure 
SO2 flows in the shell Bide and water in the tube side. The 
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as a function of liquid flow rate for module 3. 
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Figure 10. Sherwood number versus Graetz number for absorption 
of pure SO* in water. 

experimental results have been indicated for both wetted 
and nonwetted modes; the wetted mode again offers higher 
resistance to mass transfer when compared to the non- 
w ttedmode. Even though the solubility of SO2 in water 
is high and there is a dissociation reaction, the reaction 
is not fast enough to make the hquid-filled-membrane 
resistance negligible. The theoretical predictions obtained 
by numerically solving the coupled differential equations 
explain the experimental results well with a maximum 
deviation of about 13 % . The equations were also solved 
assuming that there is no reaction occurring between SO2 
andH 2 0. It was found that there was a 10% -15% increase 
in predicted SO2 concentration at the liquid outlet when 
the reaction term was included. Generally the latter 
approach describes the data much better especially in the 
wetted mode. It is worthwhile to note that a more rigorous 
analysis of SO2 dissociation in water could have been 
carried out in view of the reaction being an instantaneous 
reversible one. However, the scatter in the data (Figure 
9, nonwetted case) is considerable, and the enhancement 
due to SO2 ionization in water is known to be around 10 % 
(Danckwerts, 1970). 

Figure 10 is a pi t of Sherwood number versus Graetz 
number for the absorption of pure SO2 in water. The 
variation f Sherwood number with Graetz number is not 
very high in this case, due to the very high SO2 s lubility 
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Figure 12. Sherwood number versus Graetz number for 10% C0 2 
feed gas absorption in water. 

in water and the reaction that is occurring in the liquid 
phase. With a system having higher solubility or faster 
reaction rates, we may see a transition to gas-phase control 
(see Kaxoor (1992) for SOr-Na 2 S0 3 system). 

Figure 11 shows the results for the same configuration 
with module 2 for pure SO2 and water. At low liquid flow 
rates, when wetted and nonwetted modes of operations 
axe considered, the difference in the concentration of SO2 
at the liquid outlet is small since the gas-liquid contact 
area is high and the liquid flow rate is low so that the 
liquid gets saturated before passing the whole length of 
the module. 

So far, the cases studied utilized pure gas in the shell 
side; therefore the gas-phase resistance was negligible. 
Figure 12 shows the results for modules 2 and 3 when a 
feed gas containing 10% CO2 in N2 was flowing in the 
shell side and liquid was flowing countercurrently in the 
tube side. The results are presented as a plot of Sherwood 
number versus Graetz number. An increased resistance 
to mass transfer in the wetted mode is again apparent. 
The results of the model simulati ns show a maximum 
deviation f about 25% in the case of nonwetted mode. 
In the case f the wetted mode, the model seems to explain 
the experimental bservati ns within a maximum devi- 
ation fl8%. These deviations may have resulted from 
all f the gas n t coming into contact with the liquid in 
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Figure 14. SO> absorption in water from 1% SOy-air mixture in 
module 3. 

the modules which do not have high packing density. 
Similar trends were obtained in the case of the 40% CO«r 
balance N2 gas mixture (Figure 13); the maximum devi- 
ation in the non wetted mode was around 25% in the case 
of module 2 and around 25% in the case of module 3. In 
the wetted mode, the deviation was 10 % for module 2 and 
35% in the case of module 3. This may be attributed to 
bypassing in the shell-side flow and to the fact that the 
packing fraction is lower in the case of module 3 when 
compared to module 2. 

In these calculations the diffusion coefficient of the gas 
in the pores in the nonwetted mode was calculated from 
the transition regime of gas diffusion in porous media. 
Since rpf\ was around 0.56, the diffusion of CO2-N2 in the 
pore was considered to be in the transition regime and not 
in the Knudsen regime. It was found that if the molecular 
diffusion coefficient was used there was no difference in 
the results as CO2 absorption in water is liquid phase 
controlled. 

Figure 14 identifies the concentration of SO2 at the gas 
utlet for module 3 when 1% SO2 (balance air) mixture 
was scrubbed using water in module 3. It is seen that, as 
the liquid flow rate increases, the concentration f SO2 
decreases, as expected. When the liquid flow rate was 
increased ab ve 40 cm 3 /min, essentially n SO2 was 
obs rved in the gas utlet. For module 2 at the same gas 
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flow rate (200 cmVmin), no SO2 was observed experi- 
mentally at the gas outlet even at much lower liquid flow 
rates. This is because of the higher contacting area in the 
case of module 2 when compared with module 3. Thus it 
is observed that, with such small modules, it was possible 
to absorb all of the SO2 using simple water as the absorbent 

For the systems CO2-H2O and SO7-H2O studied, the 
rate of absorption of the gas species scrubbed using the 
same module is much higher for SO2 when compared to 
CO2. There are several reasons for this: the solubility of 
SO2 in water is 2 orders of magnitude higher than that of 
CO2 and also the dissolution of SO2 in water is accompanied 
by a reaction. Even though CO2 reacts with water, the 
rate is orders of magnitude slower when compared to SO2. 
Thus the absorption of SO2 by water is enhanced by both 
the increased solubility and the chemical reaction. 

It is useful now to identify the level of gas purification 
achieved using small hollow fibers modules having high 
surface areas. Usually most of the CO2 and SO2 were easily 
scrubbed out. In the case of SO2, removal as high as 98% 
was achieved with a liquid flow rate of 38.46 cmVmin and 
gas flow rate of 714.7 cm 3 /min in the 30-cm-long, 0.49- 
cm-diameter module 2 from a pure SO2 feed gas. When 
water is used as an absorbent, it is better to operate in a 
nonwetted mode. 

Figure 15 shows the absorption of pure SO2 and of SO2 
from SOz-air mixtures when the absorbent liquid, water, 
flows in the shell side. Pure SO2 removals as high as 99% 
were achieved with a liquid flow rate of about 40 cm 3 /min 
in a 30-cm-long, 0.45-cm-diameter module 1 with a gas 
flow rate of about 433 cm 3 / min. It is useful to know in this 
context that absorption of flue gases can be performed 
effectively using seawater (Radojevic and Tressider, 1992). 

Considerable data were taken also when either pure CO2 
or CO2-N2 mixtures flowed on the tub side while water 
flowed countercurrently n th shell side (Karoor, 1992). 
The gas abs rption rates wer always significantly 1 wer 
inmost f the systems studied c mpared to the tube-side 
liquid 11 w cases due most likely to shell-side channeling. 
A much better gas absorpti n performance with shell- 
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Figure 1 6. Comparison oiKija and HTU of a hollow fiber contactor 
to that of packed towers for CO: absorption in water, 
side liquid flow was achieved using a crossflow module 
(Karoor, 1992). These results are not being reported here. 

The pressure drop encountered in the gas-phase flow 
through the shell side was always very small for pure gas 
stream used here; mostly it was negligible. In a few cases 
the APfi encountered varied between 0.05 and 0.3 psi. When 
these gases flowed in the tube side, the AP was usually 
higher, in a few cases as much as 3-5 psi. 

The absorption performance of a hollow fiber module 
needs to be compared with that of a conventional packed 
bed absorber. To this end, several experimental obser- 
vations in a packed bed column for the CO2-H2O system 
reported in the literature are now compared to those of 
hollow fiber modules in Figure 16 using the data from the 
hollow fiber module 2 (dimensions of which are given in 
Table I). In this plot the equipment-average value of Km 
obtained from the product of experimental Koim (eq 1) 
and a is plotted as a function of superficial liquid velocity 
(a is the interfacial gas-liquid contact area per unit 
volume). To obtain the superficial liquid velocity in a 
hollow fiber module in a manner equivalent to that defined 
for a packed bed, the following definition of superficial 
velocity is used for a hollow fiber contactor: 



(30) 



60(empty cross-sectional area) 
where F\ is the liquid flow rate in cm s /min, empty cross- 
sectional area is defined as nr. 2 , and r, is the shell inner 
radius. From the figure it is observed that, as the 
superficial velocity increases, the Km also increases. Rixon 
(1948) studied the absorption of CO2 in a commercial size 
plant where the column inner diameter was 27 in. and the 
packed height was 18 ft. The packing used was 1 Vrin. 
X 1 1 / 2 -in. stoneware. The Km values obtained by them 
are quite low. Danckwerts and S harm a (1966) studied 
the absorption of CO2 in water in a column of 18-in. 
diameter and the packing material 1 V2-in. ceramic raschig 
rings. The values of Km obtained by them were also low 
when compared to those in the hollow fiber module. In 
their case, the liquid film coefficients were very high but 
the area/ volume was only 1.3 cm 2 / cm 3 , which resulted in 
low Km values. Sherwood and Holloway (1940) used their 
O2 desorption data and came up with general correlations 
f r calculating Km and HTU for packed beds with vari us 
packing materials: 
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Figure 17. Comparison of Kyi and HTU of a hollow fiber contactor 
to that of packed towers for absorption of SO2 in water. 
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where Km is the mass-transfer coefficient for liquid phase, 
h" 1 ; if l is the height of an individual liquid-film transfer 
unit, ft; L is the liquid mass velocity, lb/(ft) 2 (h); ml is the 
viscosity of the liquid, lb mass/ (ft) (h) ; pl is the density f 
the liquid, lb/ft 3 ; Dl is the diffusion coefficient for liquid 
phase, ft 2 /h (approximately 6.6 X IQ-* ft 2 /h for CO2 at 18 
°C); and a, n, and s are constants for various packings 
used. The above equations provide general estimates of 
Km and HTU values at low liquid flow rates. The above 
equations can also be used for absorption of CO2 as the 
gas-film resistance is generally neglected. As seen in Figure 
16, the Km values from such an equation are always smaller 
than those obtained in the hollow fiber contactor. Also 
most of the data for packed beds are available for low 
liquid flow rates; at higher liquid flow rates there may be 
flooding or loading. Thus, even with a simple absorbent 
liquid used in these experiments, a remarkably higher 
amount of removal was achieved with relatively small 
modules and flow velocities. The potential for increased 
removal is enormous with the application of reactive 
absorbents and appropriate operating conditions. 

Figure 16 showB also the HTU variation as a function 
of the superficial liquid velocity. The HTU in the case of 
a hollow fiber is given as 



HTU- 



(33) 



In the figure the data obtained by Cooper et al (1941) 
in a 2-in. X 2 Vw-in. steel raschig ring column show the 
highest HTU for a given superficial velocity. The HTU 
given by eq 32 (Sherwood and Holloway, 1940) is also 
plotted. As seen in the figure, the equation still yields 
higher HTU when compared to that for a hollow fiber 
contactor. The difference is especially large at lower 
superficial velocities. 

Figure 17 compares the Km values for hollow fiber 
module 3 to those of packed towers for SO2 absorption in 
water. The Km is plotted as a function of the superficial 
velocity. The Km increases as the superficial velocity 
increases. Whitney and Vivian (1949) studied the ab- 
sorption f SO2 in water in a 8-in.-diameter tower using 
1-in. ceramic packing; as bserved from the figure the Km 
btained by them were 1 w. At a superficial velocity of 
1.5 cm/a, the Km obtained in a hollow fiber module is 



around 10 times higher than that obtained in a packed 
tower. Figure 17 compares also the HTU of hollow fiber 
module 2 to that of packed towers. The data obtained by 
Whitney and Vivian (1949) and Pearson et aL (1951) are 
plotted as a function of the superficial velocity. As 
observed from the figure, very low HTU values were 
obtain d when a boll w fiber module was used f r gas- 
liquid contacting when compared to th packed towers. 

Conclusions 

For CO2 absorption in water and SO2 absorption in water 
in hydrophobic microporous hollow fiber modules, both 
liquid and gas flow rates were controlled independently, 
and there were no operational problems in either wetted 
or nonwetted modes of operation. For CO2 absorption in 
water, the wetted mode of operation offers considerably 
higher resistance to mass transfer when compared to the 
nonwetted mode of operation. It is a liquid-phase- 
controlled process; thus filling the membrane pores with 
water increases the resistance to mass transfer. Increasing 
the contact area in a given module or the packing fraction 
for a given fiber size increases the C0 2 removal capacity. 
The contact area in a hollow fiber device can be increased 
by increasing the number of fibers without any substantial 
increase in the physical size of the contacting device. When 
the absorbent liquid, water, flows in the tube side, 
simulation results from a model based on first principles 
are in good agreement with experimental observations. 
The Leveque solution also agrees well with the experi- 
mental results. The experimentally obtained K\ja values 
for the contactors are considerably larger than those for 
packed towers. Even though the mass-transfer coefficients 
obtained for hollow fiber modules in this study are similar 
to or lower than the values reported for packed beds in the 
literature, the contact area per unit equipment volume in 
the case of hollow fiber devices is many times higher than 
that for packed towers. Therefore, at a superficial liquid 
velocity of 1.25 cm/s, the K\ja value was 5 times higher 
than that for a conventional packed tower. As the liquid 
flow rate increases, the difference in the Kifl values is 
decreased. At higher flow rates, there are operational 
problems in packed tower like flooding and loading; there 
is no such problem in a hollow fiber module as long as the 
correct phase pressure conditions are maintained. The 
HTU in a hollow fiber module is several times lower than 
that in a packed tower operating under similar flow 
conditions. 

For the absorption of SO2 in pure water in a microporous 
hollow fiber contactor, the following conclusions are 
relevant. The K\ja value for a hollow fiber contactor was 
about 10 times higher than that for conventional packed 
towers for a superficial liquid velocity of 1.5 cm/s. The 
HTUs obtained were also about 10 times lower than those 
obtained for a packed tower for SO2 absorption in water. 
Thus efficient contacting is achieved in the case of a hollow 
fiber when compared to conventional contactors. In the 
case of SO2 absorption in water, removal rates as high as 
99% were attained in module 1 (a Vrin. tube containing 
102 X- 10 fibers, 30 cm long) . Thus, even in a small module 
high removal rates were attained, proving the efficiency 
of the contactor. Incorporation of SO2 dissociation reaction 
in pure water predicts 10 % -15 % higher SO2 concentration 
in the liquid outlet. The wetted mode of operation with 
pure water as absorbent offers somewhat higher resistance 
when compared with the nonwetted mode of operation. 
Thus it is advantageous to perate the module in the 
nonwetted mode. Higher contact areas per unit module 
v lume provide higher SO2 rem val efficiencies. The 
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model simulations are in good agreement with the ex- 
perimental observations for cases where the absorbent 
liquid fl wed in the tube side. 
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Nomenclature 

a » gas-liquid interfacial area/contactor volume (cm* 1 ) 
At a mass-transfer area based on surface area of gas-liquid 
contact (cm 3 ) 

Cijk m concentration of species i in the phase ; 0 * g, for gas, 
and; - 1, for liquid) and in region k (k = m, for fiber wall 
pore, k = s, for shell side, and k ■ t, for fiber bore) (mol/ 
cm 3 ) 

Caw ■ concentration of species A in liquid at the gas-liquid 
interface at the membrane wall (mol/cm 3 ) 

Cbw ■ concentration of species B in liquid at the gas-liquid 
interface at the membrane wall (mol/cm 3 ) 

C g fa, ™ concentration of the absorbed gas species at the gas 
stream inlet (mol/cm 3 ) 

Cf out ■ concentration of the absorbed gas species at the gas 
stream outlet (mol/cm 3 ) 

C ig ■ concentration of species i in the gas phase at the gas- 
liquid interface (mol/cm 3 ) 

On ■ concentration of species i in the liquid phase at the 
gas-liquid interface (mol/cm 3 ) 

C\ out ■ concentration of CO2 or SO2 in liquid outlet (mol/ 
cm 3 ) 

ACim ■ logarithmic mean concentration difference of CO2 or 
SO2 (mol/cm 3 ) 

Obi * diffusion coefficient of species B in the liquid phase 
(cm 2 /s) 

Dik * diffusion coefficient of species i in region k(k-m, for 
fiber wall pore, k ■ s, for shell side, and k • t, for tube side) 
(cmVsec) 

F t ^ ■ gas flow rate at module inlet (cmVmin) 

*V out m gas flow rate at module outlet (cm 3 /min) 

H « dimensionlesB Henry's constant as defined by eq 3 

ki » forward rate constant for S0 2 ionization (ar l ) 

Kz = equilibrium constant of a chemical reaction 

Xoim, Kl » overall mass-transfer coefficient (cm/s) 

L « effective length of the module, cm 

N( « total number of fibers in the hollow fiber module 

JVct ■ Graetz number, *rJfR>NsJ2L 

Nju - Reynolds number, 2r t i>rtp/M 

JVsc = Schmidt number, n/pD At 

Nsh ■ Sherwood number, iCoim(2rt)/D At 

Pgt Pgu m gas-phase pressure (psig) 

P«q ■ liquid-phase pressure (psig) 

r m radial coordinate position (cm) 

r, • free surface radius (see eq 6) (cm) 

r, ■ molar production rate of species i by reaction per unit 

volume (mol/(cm 3 *8)) 
r 0 , r t « outer radius and inner radius of the fiber (cm) 
r p , r» - pore diameter, shell radius (cm) 
R f u m gas constant (atm«cm 3 /(mol K)) 
Tr ■ room temperature (°C) 
v„ t v t% ■ shell aide and tube side velocity (cm/s) 
Vi « liquid flow rate (cm 3 /min) 

Greek Letters 

t - packing fraction 

My - viscosity of the fluid in the sheU side (gm/(cm-s)) 
X « mean free path (cm) 
p » density (g/cm 3 ) 



684 Ind. Eng. Chem. Res., Vol. 32, No. 4 r 1993 

r «= tortuosity of the membrane pore 
Subscripts 

A, B = species C0 2 and HS0 9 -, respectively 
i a component involved 
; = phase involved 0' • g> gas, or ; = 1, liquid) 
k * location of the Bpecies (k » m, for fiber wall pore, k * a, 
for shell side, and k » t for tube side) 
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